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HYDRODYNAMIC  RESPONSE  OF  A  RADIATING  AS  PLASMA 
EXCITED  BY  A  PROTON  BEAM 


I.  INTRODUCTION 


Recent  advances  in  pulsed  power  and  diode  technologies  have  made  it 
possible  to  deliver  intense  proton  beams  onto  a  planar  slab  with  high  effi¬ 
ciency.  For  example,  focussed  proton  beams  with  energy  of  the  order  of  an 
MeV  and  with  current  densities  of  several  kA/cm^  have  been  reported. 
Thus,  the  interaction  of  energetic  charged  particle  beams  with  planar  slab 
materials  has  become  an  active  field  of  research.  One  significant 
application  is  the  use  of  ion  beams  to  generate  x-rays  for  materials 
research,  lethality  and  vulnerability  studies. 

The  investigation  of  charged  particle  beam-matter  interaction  involves 
consideration  of  several  phenomena.  First,  there  is  the  energy  deposition 
within  the  slab  by  the  beam.  Then,  there  is  the  subsequent  hydrodynamic 
motion  of  the  slab  as  it  responds  to  this  energy  influx;  a  blow-off  region 
develops  on  the  front  (beam)  side  of  the  plasma,  while  the  back  surface 
recoils  away  from  the  beam.  As  the  slab  plasma  beats  up,  the  plasma  ions 
emit  radiation  which  can  affect  the  energy  balance  of  the  plasma.  Thus  at 
any  instant  in  the  plasma,  there  is  the  interplay  of  deposited  beam  energy, 
thermal  conduction,  hydrodynamic  motion,  radiation  transport,  plasma 
internal  energy,  magnetic  fields,  etc.  to  be  taken  into  account. 

With  the  exception  of  magnetic  and  electric  field  effects,  this  paper 
represents  an  attempt  to  treat  all  these  effects  self-consistently  In  a 
one-dimensional  simulation  of  a  proton  beam  incident  on  an  initially  solid 
density  ML  slab. 

II.  THEORETICAL  MODEL 


The  interaction  of  the  ion  beam  with  the  slab  plasma  is  strongly 
dependent  on  the  local  temperature  and  degree  of  ionization  in  the  slab. 
In  addition,  optical  pumping  and  energy  transport  by  photons  of  all 
frequencies  can  influence  the  plasma  temperature  and  degree  of  ionization, 
and  can  modify  the  hydrodynamic  response  of  the  slab.  Thus,  the  deposition 
of  the  beam,  the  hydrodynamic  evolution  and  atomic  physics  of  the  slab,  as 
well  as  the  transport  of  radiation,  must  be  calculated  self-consistently. 
Manuscript  approved  September  17,  1984. 
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Discussion  of  the  theoretical  model  can  be  separated  for  convenience 
as  follows:  (a)  hydrodynamics  and  thermal  conduction  (b)  ionization  and 
atomic  physics  (c)  radiation  emission  and  transport  and  (d)  beam  dynamics 
and  energy  deposition. 

( A)  Hydrodynamics  and  Thermal  Conduction 

The  basic  hydrodynamic  variables  of  mass,  momentum,  and  total  energy 
are  transported  in  one  dimension  using  a  numerical  scheme  with  a  sliding- 
zone  version  of  flux-corrected  transport.^  A  special  gridding  algorithm  is 
used  which  moves  zones  in  a  Lagrangian  fashion  and  adjusts  the  mesh  in 
order  to  resolve  steep  gradients  in  the  flow.  The  hydrodynamic  equations 
solved  are 


Dp 
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Dt  =  “  dx  ’ 


-  (up)  +  !rad 


+  e 


dep 


+  £-  t  n"feT  ( 


3x 


(2) 


(3) 


where  p  is  mass  density,  u  is  velocity,  P  is  pressure,  e^,  is  total  energy 
density,  efad  is  the  rate  of  energy  loss  or  gain  due  to  radiation,  Is 

the  rate  of  energy  gain  due  to  the  beam  deposition  process,  n  is  the 
thermal  conductivity,  and  V  is  the  ion  density.  The  thermal  conduction  is 
calculated  implicitly,  using  an  iterative  Crank-Nicholson  scheme. 

Since  densities  did  not  much  exceed  solid  density  in  this  study,  a 
simple  equation  of  state  was  assumed  with 

P  *  y  (eT  -  y  pu2  -  ex)  ,  (4) 


where  is  the  ion  potential  energy  due  to  ionization  and  excitation.  A 
single  temperature  model  was  employed. 
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where  m-.  is  ion  mass,  and  T  is  temperature.  The  ionization  energy,  ,  and 
effective  charge,  Z  are  calculated  from  the  ionization-radiation  equations 
which  are  explained  below.  A  single  temperature  assumption  is  valid  in  the 
solid  dense  material  as  well  as  in  the  beam  deposition  region,  since 
thermal  equilibration  times  are  short  compared  with  the  time  scales  of  the 
energy  input  and  changes  in  the  hydrodynamic  variables. 


The  local  rate  of  change  of  energy  due  to  radiation  transport  erad» 
and  that  due  to  the  beam  deposition  ,  will  be  discussed  in  Sections  C 
and  D. 


(3)  Ionization  and  Atomic  Physics 


The  ionic  populations  in  the  plasma  may  be  characterized  by  a  set  of 
atomic  rate  equations  of  the  form 


dt 


-  zj  V  f j  -  zi  wij  fi 


(6) 


where  fj  is  the  fractional  population  of  atomic  level  i,  and  Wj^  is  the  net 
reaction  rate  describing  the  transition  from  initial  state  j  to  final  state 
i.  An  equation  of  this  type  is  constructed  for  each  of  the  atomic  levels 
included  in  the  model. 

For  sufficiently  dense  plasmas  of  the  sort  we  expect  to  model  with  the 
ion-beam  code,  the  effective  populating  and  depopulating  rates  are 
generally  fast  compared  with  the  hydrodynamic  response.  An  equilibrium 
assumption  can  be  justified,  which  involves  dropping  the  explicit  time 
dependence  in  equation  (6).  The  plasma  Is  then  said  to  be  in  collisional- 
radiative  equilibrium  (CRE),^  whereby  the  plasma  ionization  state  responds 
instantaneously  to  changes  In  hydrodynamic  quantities. 

The  rate  coefficients  that  are  used  to  calculate  the  populating  and 
depopulating  rates,  Wji»  are  calculated  using  various  scattering 
techniques.  The  processes  included  In  this  calculation  and  the  methods 
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used  in  calculating  the  corresponding  rate  coefficients  are  summarized 

below. 

(1)  Collisional  ionization  -  exchange  classical  impact-parameter  (ECI?) 
methods'^  (the  effect  of  autoionizing  resonances  on  the  cross  sections  has 
been  ignored). 

(2)  Photo ionization  -  hydrogenic  approximation  with  Karzas-Latter-Gaunt 
factors^*®. 

(3)  Dielectronic  recombination  -  the  detailed  calculations  of  Jacobs  et 
o 

al.  are  used. 

(4)  Collisional  excitation  -  Coulomb-Born  distorted-wave  approximation 
including  exchange  effects*®,  or  the  semiclassical  impact-parameter  (SCI) 
technique**. 

(5)  Spontaneous  radiative  decay  -  oscillator  strengths  are  taken  from 
several  calculations  and  measurements *^. 

(6)  Photoexcitation  -  oscillator  strengths  used  are  those  quoted  above  to 
determine  optical  depths  (see  the  next  section  on  radiation  transport). 

Finally,  collisional  and  radiative  recombination,  collisional  deexcitation, 
and  stimulated  emission  are  all  calculated  as  the  detailed  balance  of  the 
corresponding  opposite  rate  listed  above. 

Once  the  set  of  rate  equations  (including  the  radiation  transport)  has 
been  solved  for  the  level  populations,  f^,  the  electron  density  can  be 
calculated , 
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E  *  f^ 
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N. 
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where  z^  is  the  ionic  charge  of  level  i  and  Nj  is  the  total  ion  density. 


The  ionization  and  excitation  energy  can  also  be  calculated  by 


"i 


fi  V 


(8) 


where  is  the  energy  of  Level  i,  measured  from  the  ground  state  of  the 
neutral  atom. 


(C)  Radiation  Emission  and  Transport 

Radiation  emission  from  a  plasma  and  its  opacity  are  dependent  on  the 
local  atomic  level  population  densities.  Except  for  optically  thin 
plasmas,  however,  the  level  populations  depend  on  the  radiation  field, 
since  optical  pumping  via  photoionization  and  photoexcitation  can  produce 
significant  population  redistribution.  Thus,  the  ionization  and  radiation 
transport  processes  are  strongly  coupled  and  must  be  solved  self- 
consistently.  In  this  model,  an  iterative  procedure  J  is  used,  where  level 
populations  are  calculated  using  the  radiation  field  from  the  previous 
iteration,  then  using  these  populations  to  calculate  a  new  radiation  field 
and  recalculating  populations  until  convergence  is  reached. 

Three  distinct  radiation  transport  schemes  have  been  developed,  and 
can  be  used  interchangeably  in  the  code:  A  probabilistic  scheme,  a 
multifrequency  scheme  and  a  hybrid  method.  Descriptions  of  these  transport 
methods  have  appeared  elsewhere;  in  this  report,  their  general  properties 
will  be  outlined. 

The  probabilistic  model ^  forms  local  angle  and  frequency  averaged 
escape  probabilites  for  each  emission  line  and  for  each  bound-free 
process.  Free-free  radiation  is  treated  with  a  multifrequency  formalism. 
The  radiation  transport  and  emission  spectra  are  calculated  from  these 
escape  probabilities.  The  method  is  cost-effective,  can  treat 
comprehensive  atomic  models  and  provides  good  overall  energetics,  but 
cannot  calculate  certain  spectral  details  and  breaks  down  at  very  high 
densities. 

The  multi frequency  model ^  solves  the  equation  of  radiation  transport 
at  a  large  number  of  discrete  frequencies,  providing  resolution  of  emission 
lines,  recombination  edges  and  absorption  edges.  It  provides  accurate 
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radiation  transport  at  high  density,  and  gives  spectral  details  such  as 
self-absorption  features.  However,  a  large  number  of  frequencies  is 
required  to  provide  adequate  resolution,  and  it  tends  to  be  a  relatively 
costly  method. 

The  hybrid  model  uses  the  multifrequency  formalism  to  transport  the 
continuum  (bound-free  and  free-free)  radiation  and  frequency-integrated 
escape  probabilities  to  transport  the  lines.  Continuum  opacities  are 
interpolated  from  the  multifrequency  mesh  and  folded  into  the  line 
transport  calculations.  The  line  opacities  are  assumed  to  have  negligible 
Impact  on  the  continuum  energetics. 

In  all  of  the  models  outlined  above,  local  ionization  state-dependent 
inner-shell  opacities  are  included,  since  these  processes  are  very 
important  in  the  cool,  dense  plasma  regions.  Inner-shell  photoionization 
cross  sections  for  the  neutral  element  are  taken  from  the  fits  by  Biggs  and 
Lighthill,^  and  the  positions  of  the  ionization-dependent  absorption  edges 
are  taken  from  the  Hartree-Fock  calculations  of  Clementi  and  Roetti. ^ 

The  local  rate  of  energy  change  in  zone  j,  due  to  radiation  transport 
is  given  by 

b  - '  l  <fpj  - 1  Skj  FPk>  (9> 

where  Fp^  is  the  rate  of  energy  loss  in  zone  k  due  to  a  discrete  radiative 
process  (or  frequency  group)  P,  and  Cp^j  *s  th®  radiative  coupling  of  zone 
k  to  zone  j  for  that  process.  The  couplings  are  functions  of  opacity, 
integated  over  process  and  photon  path.  In  the  probabilistic  model,  a 
matrix  of  couplings  must  be  computed  for  each  bound-bound,  bound-free  and 
free-free  process;  for  the  multifrequency  model,  they  must  be  computed  for 
each  discrete  frequency.  In  this  way,  the  net  cooling  and  heating  by 
radiation  emission  and  absorption  among  the  various  zones  of  the  plasma  is 
accounted  for  accurately. 

(0)  Beam  Dynamics  and  Energy  Deposition 


We  assume  that  on  the  timescale  of  the  slab  response  (several 
nanoseconds),  the  beam  maintains  uniform  flow,  that  is 
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(10) 


where  and  Ug  are  the  local  mass  density  and  axial  velocity  of  the 
beam.  The  mass  flux  m  is  assumed  constant  up  to  the  point  the  beam  is 
stopped.  J  is  the  beam  current  per  unit  area;  mg  and  Zg  are  the  mass  and 
charge  of  a  beam  ion. 


For  a  relativistic  beam  with  negligible  thermal  spread,  the  kinetic 
energy  of  a  beam  ion  is  related  to  the  local  velocity  by 
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The  energy  is  also  related  to  the  stopping  power  by 
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where  EQ  is  the  initial  energy,  x"  is  the  distance  into  the  slab,  and  the 
total  stopping  power  is  given  by  the  sum  of  the  ion-free  electron,  ion- 
bound  electron  and  ion-plasma  ion  interactions. 
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1.  Interaction  With  Bound  Electrons 


Two  widely  used  models  for  calculating  stopping  cross  sections  are  the 
LSS  model,  which  is  valid  at  low  velocity,  and  the  Bethe  theory,  which  is 

valid  at  high  velocity.  These  are  complementary  models  and  can  be  used  to 

1  ft 

span  the  entire  projectile  velocity  range  [e.g.,  Mehlhorn1-  ]  . 

The  chief  difficulty  in  applying  the  Bethe  theory  is  obtaining  the 
mean  excitation  energy  T.  The  calculation  of  this  quantity  is  a  tedious 
exercise,  and  various  scaling  schemes  have  been  proposed  to  estimate  it 
[Mosher^,  Mehlhorn^®) .  Also,  alternate  methods  of  calculating  bound 
electron  stopping  power  have  been  developed  which  do  not  involve  this 
quantity.  One  of  these  is  the  local  oscillator  model  [Rogerson  ,  Nesbet 
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and  Ziegler  ] .  This  model  involves  an  integral  over  the  atomic  (or  ionic) 
electron  density  and  thus  depends  on  the  accuracy  of  the  electron  distri¬ 
bution  used  in  the  calculations.  For  cold  targets,  a  Thoraas-Fermi  density 

has  been  found  to  yield  accurate  stopping  powers  for  cold  target  ions 
no 

■ Rogersofl"  . . 

-)  o 

In  the  local  oscillator  model  (LOM)  representation  of  the  stopping 
cross  section,  the  spatial  change  of  energy  given  by 

(f)b  -  (POO/mj)  Sb(V)  (14) 


where  Sb(V)  is  the  stopping  cross  section  in  erg-cm  . 

2  4  r 

4ttZ  e  f  R  ,  ^ 

Sb(V)  *  - ^ -  J  4tt  r“  p  (r)  Kq(t)  dr,  (15) 

mV"  o 

where  Zg  is  the  effective  charge  of  the  projectile  ion,  R  is  the  ion 

radius,  e  and  m  are  the  electron  charge  and  mass,  V  is  the  projectile 

velocity,  p(r)  is  the  local  bound  electron  density  in  the  target  atom,  and 

X  (t)  is  a  modified  zeroth  order  Bessel  function, 
o 

t  *  -ft  (d  (r)/mV^  ,  (16) 

o 

where  "ff  is  Planck's  constant  divided  by  2it  ,  and  ^^(r)  is  the  local  plasma 
frequency  at  radius  r  within  the  atom. 

There  are  three  assumptions  underlying  the  LOM^l.  The  first  is  that  a 
loss  function  can  be  defined,  dependent  only  on  the  local  electron  density 
in  the  target  atom.  The  second  is  that  the  longitudinal  dielectric 
response  can  be  represented  by  e(ui)  with  a  single  zero  at  w  -  subject 

to  the  high  frequency  condition 

<•  (u>)  '  1  -  w^/a>2  (17) 

appropriate  to  free  electrons.  The  third  assumption  is  that  the  induced 
polarization  charge  is  spread  out  from  the  ion  trajectory  to  some  finite 
radius  of  the  order  of  the  de  Broglie  wavelength  -ff/mV.  This  last 
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assumption  is  justified  by  the  adiabatic  argument  that  energy  will  be 
transferred  only  to  electrons  with  velocity  less  than  V.  A  wave  packet 
with  momentum  of  order  mV  would  have  a  spatial  spread  of  order  di/mV. 

The  electron  density  o(r)  is  taken  from  the  'Hiomas-- ermi  ( model 
proposed  by  Zink.~~ 

In  a  plasma  target,  the  bound  atomic  electrons  can  be  screened  from 
the  projectile  ions  by  the  free  electrons.  Thus,  the  LOM  must  be  modified 
to  take  this  shielding  effect  into  account.  For  an  electron  in  an  isolated 

atom,  the  maximum  impact  parameter  is  V/w,  where  to  is  a  characteristic 

2  2 

frequency  of  motion  .  Plasma  screening  limits  this  parameter  to  the  Debye 
length  D.  The  expression  for  x  can  be  rewritten 


r  =  (VmV)  Ooo(r)/V).  (18) 

'.whenever  V/u)  (r)  >  D,  t  is  taken  to  be 

i  =>  -fi/mVD 

Thus,  plasma  shielding  of  the  bound  electrons  is  taken  into  account  by 
limiting  the  argument  of  the  Bessel  function. 

Some  criticism  of  the  use  of  free-electron-gas  approximations  for 
bound  electron  stopping  power  (which  includes  the  LOM)  have  been 
expressed"^;  it  is  argued  that  the  use  of  a  free  electron  gas  model  to 
calculate  the  Bethe  mean  excitation  energy  parameter  I(Z,q)  -  where  Z  is 
the  nuclear  charge,  and  q  is  the  ionic  charge-yields  the  wrong  Z-scaling 
for  hydrogenic  ions;  i.e.,  this  model  yields 

I(Z,Z-1)  =■  Z3/2  1(1,0) 


instead  of  the  correct 


I(Z,Z-1)  -  Z2  1(1,0). 

Comparison  of  I(Z,q)  for  A1  derived  from  the  free-electron-gas  approxi¬ 
mation  with  I(Z,q)  derived  from  detailed  calculation  of  A2.  electron  power 
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cor  protons  show  an  increasing  divergence  with  increasing  q,  with  the 
f ree-electron-gas  I(Z,q)  being  lower. 

Recently,  McGuire,  Peek,  and  Pitchford"^  (MPP)  published  detailed 

calculations  of  stopping  power  cross  sections  of  protons  by  A2  ions.  By 

using  the  generalized-oscillator- strength  formulation  of  the  Born 
approximation,  they  calculated  tables  of  stopping  cross  sections  for 
protons  interacting  with  A£  ions  of  charge  q,  where  0  <  q  <  11,  in  the 

energy  range  0. 1  MeV  to  100  MeV. 

The  MPP  cross  sections  for  neutral  Afc  agreed  very  well  with  local 

20 

oscillator  model  calculations  ,  but  for  ionized  M,  the  MPP  cross 
sections  were  generally  lower,  (especially  at  lower  energies)  than  the 

local-oscillator-model  predictions.  Figures  1  and  2  show  comparisons  of 
these  cross  sections  for  Al+*  and  Al+^;  these  figures  indicate  generally 
the  relative  comparisons  of  the  two  sets  of  cross  sections. 

By  use  of  least- squares  methods,  these  MPP  data  were  fitted  by  a 
function  of  the  form 


Fq(y)  3  exp  (Aq  +  Bqy  +  Cqy2  +  Dqy3)  , 


where  y  =  In  (E/EQ),  and  EQ  =  0.1  MeV.  F  (y)  gives  stopping  cross  sections 
—  IS  7  ^ 

in  units  of  10  eV  cm  /atom.  Table  1  gives  the  fitting  parameters  Aq, 
3q,  0^  and  Dq  for  each  value  of  q.  Table  2  gives  the  MPP  data  and  the 
results  of  Eq.  (20).  All  the  fits  were  within  7%,  except  for  q=3,  where 
the  maximum  deviation  was  9%.  Thus,  Eq.  (20)  yielded  good  fits  to  the 
calculated  cross  section. 


In  the  treatment  of  the  beam- target  interaction,  the  average  ionic 
charge  Z  in  a  given  cell  is  used  to  obtain  stopping  power.  For  non-integer 
values  of  Z,  Eq.  (20)  is  used  for  q  and  q+1,  where  q  is  the  integer  part 
of  Z,  and  the  stopping  cross  section  for  Z  is  obtained  by  linear 
interpolation  between  Fq(y)  and  Fq+^(y). 

For  q=0,  q=l  and  q=2,  the  MPP  data,  and  hence  the  fits,  are 
monotonically  decreasing  for  energies  _>.  0*  1  MeV.  For  these  cases,  a  local- 
oscillator  model  (L0M)"  is  used  to  calculate  bound  electron  stopping  power 
for  energies  below  0.1  MeV.  These  L0M  cross  sections  are  adjusted  to  match 


$ 


the  MP?  data  at  0.  1  MeV,  and  the  cross  sections  at  lower  energy  are  scaled 

accordingly.  The  MPP  data  for  q  ■  8,  9,  and  10  are  also  monotically 

decreasing  for  E  _>  0.  1  MeV,  but  at  these  charge  levels,  the  stopping  power 

is  dominated  by  the  f ree-electron  contribution;  hence,  no  adjustments  to 

*  ■ 

the  bound  electron  stopping  cower  at  energies  less  than  0.  1  MeV  are  made. 
The  fits  are  just  extended  below  0.1  MeV. 

7or  the  other  q  values  between  3  and  11,  the  MPP  data  peak  in  the  0.1 
MeV  -  100  MeV  range;  hence  the  fits  are  extended  below  0.1  MeV  to  obtain 
low  energy  cross  section  data.  For  q  >  11,  the  bound  electron  stopping 
power  is  set  to  some  low  value,  since  free  electron  stopping  dominates. 

From  an  energetics  point  of  view,  errors  in  the  low  energy  regime  are 
not  as  critical  as  at  higher  energies.  For  example,  a  proton  with  an 
initial  energy  of  one  MeV  has  lost  90%  of  its  energy  and  is  near  the  end  of 
its  range  by  the  time  it  degrades  to  a  0.1  MeV  proton.  Therefore,  errors 
in  stopping  power  below  0.  1  MeV  will  not  affect  energy  deposition  profiles 
very  much. 

(2)  Interaction  With  Free  Electrons 


For  heated  target  materials,  the  atoms  become  ionized,  and  stopping 
due  to  free  electrons  must  be  considered.  The  free  electron  stopping  power 
is  calculated  from 


(21) 


with 


Sf(V) 


2Tte4  Z  2  Z 


F(C>  {  In  (1  +  2j)  +  in  (1  +  40  } 


(22) 


where  Z  is  the  average  charge  of  the  target  ions,  D  is  the  Debye  shielding 
length,  and 


5  -  (mBV2/2kT)1/2  , 
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(23) 


where  kT  is  the  electron  temperature,  and 


F(  ?)  =  erf  (•)  -  2  ;  e  /  /t  . 


The  quantity  b  is  the  minimum  impact  parameter  for  electron-ion 


scattering  and  is  given  by 


H _  L 

»  J  » 


i.e. ,  the  maximum  of  either  the  classical  or  quantum-mechanical  impact 
parameter  defined  by  the  uncertainty  principle. 

The  first  term  in  S^(V)  is  the  short  range  ion-electron  binary- 

0  f. 

encounter  scattering  term  and  is  taken  from  the  work  of  Campbell,'-  who 
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adapted  it  from  Brueckner  and  Brysk.  The  second  term  is  the  polarization 

OQ 

term  and  is  taken  from  Pines  and  3ohm.  For  distances  larger  than  D,  the 
plasma  acts  as  a  continuous  medium,  and  distant  collisions  cause  loss  of 
energy  by  the  excitation  of  plasma  oscillations,  which  appear  as  an 
oscillating  wake  behind  the  projectile. 


(3)  Interaction  with  Plasma  Ions 


For  stopping  power  due  to  the  plasma  ions,  a,  expression  taken  from 
Campbell^  and  Mehlhorn^  is  used; 


'S-'i "  ■mr  si<v)> 


where 


where 


U  2-2 
4me  Z* 

S1(V) - \ -  in  A  F(0 

m^V 


?2  -  mIV2/2kTi  , 


4 - - 


e 


in  A  =  ( 1  ' 2 )  in  (I  +  b  “  /  b  7  ) 

max  min 


b  is  set  equal  to  D,  and  b  ..  is  given  by 


Z  Ze2 

b»i„-MAXl-5^--5snrl 

umpV  M  P 

where  nip  is  the  proton  mass,  and 


*3  mT 


(4)  Enhanced  Stopping  Power  From  Collective  Effects 


The  stopping  power  of  a  dense  ion  beam  can  differ  from  the  sum  of 
single  particle  stopping  powers.  One  way  in  which  this  can  occur  is 
through  phase  mixing  of  the  polarization  wakes  produced  in  the  ambient 
medium  by  the  beam  particles.  The  magnitude  of  this  particular  effect  is 
similar  to  that  given  by  McCorkle  and  Iafrate  for  electron  beams. 


Sb  “  Sq  ( 1  +  2tt  rig  aJ/3) 


where  is  the  enhanced  stopping  power,  SQ  is  the  single  particle  stopping 
power,  nB  is  the  beam  number  density  and  a  is  the  ratio  of  the  beam 
velocity  and  local  plasma  frequency 
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For  the  case  considered  below,  a  1  MeV  proton  beam  with  a  flux  n^  v^  “  10 

—  0  —  1 

cm  “  sec  on  a  solid  density  aluminum  target,  the  correction  to  the 
stopping  power  is  quite  small. 


III.  RESULTS 


Two  treatments  of  the  proton  beam-  A l  target  interaction  have  been 
performed.  The  first  used  the  LOM  to  calculate  bound  electron  stopping, 

whereas  the  second  calculation  utilized  the  MPP  stopping  cross  sections. 

Figures  1  and  2  indicate  that  these  two  sources  give  different  stopping 

powers,  particularly  at  lower  proton  energies.  Thus,  two  simulations  were 
done  to  test  the  effects  of  these  different  cross  sections.  On  the 
assumption  that  the  MPP  data  are  more  accurate,  the  emphasis  will  be  on  the 
results  from  this  treatment.  Some  comparisons  of  results  from  the  two 
simulations  will  be  given. 

The  interaction  of  a  monoenergetic  1  MeV  proton  beam  with  a  planar 
aluminum  slab  of  15  pm  thickness  is  treated  in  this  investigation.  The 
beam  is  assumed  to  consist  of  a  square-shaped  pulse  of  10  nanosecond 

duration  with  a  flux  of  10^  protons /cm^-sec.  The  beam  intensity  on  target 
is  1.6x10^  W/cm^,  which  is  comparable  to  the  intensities  available  with 
current  devices. 

A  slab  thickness  of  15  pm  is  chosen  because  this  thickness  corresponds 
to  the  range  of  1  MeV  protons  in  cold  solid  density  A1  .  Hence,  the  beam 
is  initially  totally  stopped  within  the  target  and  deposits  most  of  its 
energy  near  the  back  surface  due  to  the  Bragg  peak  in  the  stopping  cross 
sections.  However,  as  may  be  seen  from  Figs.  1  and  2  and  Table  2,  the 
bound  electron  stopping  power  decreases  as  the  A1  ion  charge  increases; 
hence,  as  the  target  temperature  rises,  and  the  ion  charge  increases,  the 
beam  penetrates  completely  through  the  target  until  the  free  electron 
population  rises  sufficiently  to  stop  the  beam  entirely  within  the  plasma. 

A  point  to  be  noted  also  in  this  regard  is  that,  for  consistency,  the 
average  charge  obtained  from  solving  the  CRE  equation  is  used  throughout 
the  calculation.  Continuum  lowering  is  not  taken  into  account.  This  means 
that,  near  solid  density,  the  average  charge  Z,  and  hence  the  electron 
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density,  is  probably  low.  The  free  electron  contribution  to  stopping  power 

is  somewhat  underestimated  at  early  times  in  the  plasma  evolution. 

However,  on  a  time  scale  of  the  order  of  nanoseconds,  the  plasma  density  is 

two  orders  of  magnitude  below  solid  densitv  p  and  the  CRE  Z  is  valid. 

s 

Figures  3  and  4  show  density  and  temperature  profiles  .at  different 
times  as  the  plasma  evolves.  As  may  be  noted  from  these  figures,  the 
plasma  spreads  and  becomes  more  uniform  in  density  and  temperature  as  time 
increases. 

Figure  5  compares  density  and  temperature  profiles  at  one  nsec  from 
the  treatment  using  MPP  data  with  the  treatment  using  the  LOM  model.  Use 
of  the  smaller  MPP  cross  sections  results  in  deeper  penetration  by  the 
beam.  Hence,  the  backside  density  is  lower,  and  the  temperature  higher, 
than  that  which  results  from  use  of  the  LOM  stopping  power.  The  LOM  result 
shows  a  sharper  density  gradient  near  the  rear  surface.  Away  from  the 
backside,  the  density  and  temperature  results  are  reversed;  the  LOM  curves 
are  hotter  and  less  dense  due  to  increased  energy  deposition  in  that  part 
of  the  plasma. 

Figure  6  compares  density  and  temperature  profiles  from  the  two 
calculations  at  about  3  nsec.  The  MPP  density  curve  is  much  more  uniform 
now,  as  is  the  LOM  curve,  except  that  a  steep  density  gradient  at  the  rear 
surface  is  still  maintained.  The  LOM  still  yields  higher  temperature  and 
lower  density. 

Figure  7  compares  density  and  temperature  profiles  at  10  nsec  from  the 
two  calculations.  At  this  late  time,  the  two  profiles  are  similar.  By 
this  time,  beam  stopping  is  dominated  by  free  electron  stopping,  and  the 
difference  in  the  bound  electron  stopping  power  models  has  become 
irrelevant. 

Figures  8  through  11  show  calculated  spectra  from  the  MPP  simulation 
emitted  from  the  front  and  rear  sides  of  the  plasma.  In  Fig.  8,  at  0.5 
nsec,  the  rear  side  at  temperatures  near  40  eV  is  about  twice  as  hot  as  the 
front  surface  due  to  peaking  of  the  stopping  power  near  the  rear  surface. 
This  results  in  a  more  intense  rearside  spectrum.  As  may  be  seen  from  Fig. 
4,  this  situation  has  reversed  at  later  times,  and  at  4.31  nsec,  the  front 
side  is  emitting  a  much  more  intense  spectrum  than  the  rear  side.  The  peak 
intensity  on  the  front  side  has  risen  two  orders  of  magnitude  and  shifted 
to  higher  photon  energies.  The  rear  side  intensity  has  also  increased,  but 


less  dramatically.  Some  K  shell  emission  in  the  region  around  one  to  two 
keV  can  be  noted. 


Front  and  rearside  spectra  at  about  10  nsec  (when  the  beam  is  shut 
off)  are  shown  in  Fig.  10.  The  front  side  spectrum  is  very  similar  to  the 
4.31  nsec  spectrum,  but  the  near  side  intensity  has  increased  as  the  plasma 
continues  heating  and  becomes  more  uniform  [see  Figs.  3  and  4].  K  shell 
emission  is  very  prominent  here.  Figure  11  shows  spectra  at  18  nsec  when 
the  plasma  has  cooled  considerably  and  has  become  very  uniform.  The 
spectra  from  both  sides  are  very  similar  and  have  been  reduced  from  the  10 
nsec  levels. 

For  the  sake  of  comparisons,  front  and  rear  spectra  at  one  nsec  from 
the  MPP  simulation  and  from  the  LOM  simulation  are  shown  in  Figs.  12  and 
13.  As  in  Fig.  8  for  the  MPP  simulation  at  0.5  nsec,  Fig.  12  shows  the 
backside  emission  to  be  more  intense.  As  can  be  noted  from  Fig.  5,  the 
rear  surface  is  hotter  and  slightly  more  dense  than  the  front  surface.  In 
Fig.  13,  the  situation  is  reversed,  with  the  front  surface  being  the 
stronger  emitter  than  the  rear  surface  in  the  LOM  case.  From  Fig.  5,  for 
this  situation,  the  front  surface  is  hotter  and  much  less  dense.  In  Figs. 
12  and  13,  the  front  side  spectra  are  very  similar  due  to  similar  density 
and  temperature,  but  t.h rear  MPP  spectra  are  much  more  intense  due  to 
higher  temperature  and  lower  density. 

Figure  14  shows  front  and  rear  spectra  at  about  10  nsec  from  the  LOM 
simulation.  Comparison  with  Fig.  10  shows  these  spectra  to  be  very  similar 
to  the  MPP  spectra  at  the  same  time.  This  results  agrees  with  Fig.  7, 
where  it  was  seen  that  the  10  nsec  density  and  temperature  profiles  from 
the  two  calculations  are  very  similar.  Thus  it  may  be  concluded  that  at 
late  times,  the  differences  in  the  two  atomic  electron  stopping  powers  do 
not  contribute  substantially  to  the  results,  particularly  diagnostic  tools 
like  spectra. 

Figure  15  shows  the  energy  history  of  the  plasma  and  the  partitioning 
of  that  energy.  At  early  times,  the  radiated  power  is  an  insignificant 
fraction  of  the  beam  energy,  but  at  late  times,  it  becomes  a  significant 
portion  of  the  energy  budget.  At  about  12  nsec,  it  equals  then  surpasses 
the  total  energy  (sum  of  kinetic,  potential,  and  thermal  energies)  in  the 
plasma;  beyond  this  point,  more  of  the  deposited  energy  has  been  radiated 
away  than  remains  in  the  plasma. 


To  investigate  the  importance  of  radiation  transport,  the  MPP 
simulation  was  repeated,  but  with  radiation  transport  omitted.  Fig.  16 
shows  temperature  and  density  profiles  at  3  nsec  with  and  without 
radiation.  Including  radiation  effects  produces  a  more  unif ormly-dense 
plasma  in  this  case.  The  densities  are  not  very  different,  except  that  the 
profile  without  radiation  shows  a  dip  in  the  center  and  a  much  less  dense 
blow  off  region  in  the  front.  The  radiationless  temperature  profile  also 
attains  a  higher  maximum  value. 

Figure  17  shows  similar  comparison  at  10  nsec.  Both  density  profiles 
exhibit  gradients  near  the  rear  surface,  but  the  radiationless  density 
gradient  is  much  steeper,  and  has  a  lower  density  blow  off  region  in  the 
front.  Because  there  is  no  energy  loss  due  to  radiation,  much  higher 
temperatures  result. 

Figure  18  shows  the  temperature  and  density  of  the  rearside  of  the 
plasma  as  a  function  of  time  when  radiation  is  included  and  when  it  is 
omitted.  The  curves  are  very  similar  until  about  4  nsec;  after  this  time 
the  rear  side  is  denser  and  cooler  when  radiation  is  omitted.  Since,  at 
these  times,  the  beam  is  stopped  in  the  plasma  well  before  it  reaches  the 
back  side,  the  large  temperature  difference  between  4  nsec  and  10  nsec 
indicates  that  radiation  transport  plays  a  significant  role  in  heating  the 
rear  side.  This  agrees  with  a  conclusion  reached  in  a  study  of  the 
interaction  of  a  10iJ  W/cm  laser  pulse  at  1.06  pm  wavelength  on  8.0  pm  A* 
[Duston,  et.al^]* 

The  density  and  temperature  rise  around  9  nsec  in  the  radiationless 
case  is  due  both  to  a  shock  reaching  the  rear  surface  and  to  deeper 
penetration  by  the  beam.  This  deeper  penetration  results  from  generally 
lower  plasma  densities  at  later  times  than  when  radiation  is  included 
(e.g.,  see  Fig.  17).  Rearside  heating  at  this  time  is  due  primarily  to 
convection.  In  the  other  treatment,  the  beam  is  still  being  stopped  near 
the  center  of  the  plasma,  but  radiation  transport  is  keeping  the  rear  side 
at  a  high  temperature  even  after  the  beam  turns  off  at  10  nsec. 

Fig.  19  shows  the  time  behavior  of  the  density  and  temperature  of  the 
front  side  of  the  plasma.  Here  the  effects  of  radiation  transport  are  more 
pronounced;  the  radiationless  case  produces  much  higher  temperatures  and 
much  lower  densities.  The  marked  drop  in  temperature  after  10  nsec  is  due 
to  radiative  cooling  in  the  case  where  it  is  included. 


rig.  20  gives  Che  energy  history  of  the  plasma  along  with  the  energy 
partitioning  for  the  radiationless  case.  The  plasma  energy  is  equal  to  the 
deposited  energy  at  all  times.  The  thermal  energy  is  much  higher  than  in 
rig.  15  for  the  radiation-hydro  case;  this  indicates  that  radiation  cooling 
occurs  at  the  expense  of  the  plasma  thermal  energy.  This  is  also  indicated 
by  the  much  higher  temperatures  in  the  radiationless  curves  in  Figs.  16- 
17.  The  ionization  energies  are  not  very  different  in  the  two 
calculations.  The  radiationless  case  has  higher  kinetic  energy  at  later 
times;  this  accounts  for  the  greater  widths  of  the  hydrodynamic  plasma. 

These  comparisons  of  the  radiationless  and  hydro-radiation  transport 
cases  clearly  establish  the  importance  of  radiation  transport  effects  in 
hot  plasmas.  Any  realistic  treatment  of  hot  plasmas  must  Include  self- 
consistent  radiation  effects. 

IV.  COMPARISON  WITH  LASER-TARGET  INTERACTION 

It  is  of  interest  to  make  some  qualitative  comparisons  with  a 
simulation  of  a  laser-aluminum  target  interaction  at  roughly  the  same 
incident  energy  intensity.  Duston,  et.  al.  studied  the  interaction  of  a 
1.06  micron  laser  beam  at  an  intensity  of  10*  W/cm  and  a  long  pulse 
length  (3  nsec  full  width  at  half  maximum)  with  an  8  pm  thick  aluminum 
foil.  Fig.  21  is  taken  from  that  study. 

In  the  laser  case,  after  the  initial  blow  off  plasma  is  formed,  the 
beam  does  not  penetrate  into  the  target  beyond  the  critical  density 
surface,  whereas  the  ion  beam  penetrates  deep  into  the  target,  at  least  at 
early  times.  In  both  cases,  the  dense  part  of  the  plasma  at  later  times  is 
heated  predominantly  by  radiation  transport  rather  than  thermal 
conduction.  The  laser  case  shows  much  lower  densities  and  higher 
temperatures  on  the  front  side  of  the  plasma,  whereas  on  the  backside  the 
laser  case  shows  much  higher  density  and  lower  terapertures.  Also,  by 
comparing  Figs.  3  and  21,  it  can  be  seen  that  the  laser  beam  calculation 
produces  much  sharper  density  and  temperature  gradients.  These  results 
substantiate  the  view  that  proton  beams  are  more  efficient  than  laser  beams 
for  heating  up  the  interior  and  backside  of  targets;  this  is  due  to  deeper 
penetration  into  the  target  by  the  ion  beams. 


V.  SUMMARY 


A  fully  self-consistent,  one-dimensional  treatment  of  the  interaction 
of  a  one-MeV  proton  beam  with  an  aluminum  slab  has  been  presented.  A  fully 
self-consistent  hydrodynamic-ionization-radiation  transport  model  was 
used.  The  atomic  populations  as  functions  of  density  and  temperature  were 
obtained  from  a  collisional-radiative-equilibrium  assumption.  Radiation 
transport  was  calculated  using  a  hybrid  scheme  that  combines  a  multi¬ 
frequency  formalism  for  the  continuum  with  a  probability-of-escape  method 
for  the  lines. 

Profiles  of  the  plasma  density  and  temperature  at  various  times  in  the 
plasma  evolution  have  been  presented.  Calculated  spectral  emissions  from 
the  front  and  rear  surfaces  at  varius  times  have  also  been  presented. 

A  simulation  involving  bound  electron  stopping  power  derived  from  a 
local-oscillator-model  (LOM)  has  been  compared  with  a  simulation  involving 
the  detailed  stopping  cross  sections  of  McGuire,  Peek.,  and  Pitchford  (MPP), 
which  are  generally  lower  at  lower  proton  energies.  At  early  times  (of  the 
order  of  a  nanosecond),  the  different  stopping  powers  produce  different 
profiles  and  spectra  (see  Figs.  5-6  and  Fig.s  12-13),  but  at  later  times, 
free  electron  stopping  dominates  bound  electron  stopping  (for  this 
particular  combination  of  beam  power  and  slab  thickness  and  the  differences 
in  the  two  simulations  disappear  (see  Figs.  7,  10,  and  14). 

Comparisons  with  calculations  omitting  radiation  show  that  radiation 
transport  can  significantly  affect  plasma  evolution,  leading  to 
substantially  different  temperature  and  density  profiles,  and  must  be 
included  in  any  self-consistent  treatment  of  ion  beam-target  interactions. 

Qualitative  comparisons  of  the  ion  beam- slab  interaction  with  a 

O  1 

treatment-^  of  a  laser  beam-slab  interaction  shows  that  the  laser  produces 
higher  temperatures  and  lower  densities  on  the  front  (beam)  side  of  the 
plasma,  while  the  proton  beam  is  more  efficient  at  heating  up  the  back  side 
due  to  deeper  penetration  into  the  target. 

Since  these  treatments  of  the  beam-target  interaction  are  one¬ 
dimensional,  energy  flow  normal  to  the  direction  of  the  beam  has  been 
neglected;  in  addition,  self-field  effects  have  not  been  considered.  These 
effects  could  alter  the  range  significantly,  and  alter  the  magnitude  and 


shapes  of  the  plasma  profiles  given  here.  A  two-dimensional  radiation 
hydrodynamics  code  will  be  used  to  investigate  some  of  these  effects  in  a 


separate  paper. 
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Table  1  —  Least-square-fitting  parameters  to  the  MPP  stopping  powers* 
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*  Numbers  in  parentheses  indicates  powers  of  ten 
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Fig.  14.  Front  and  rear  spectra  at  10  nsec  (LOM  cross  sections). 
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Fig.  16.  Comparison  of  temperature  and  density  profiles  at  3  nsec 
including  radiation  transport  with  temperature  and  density 
profiles  omitting  this  effect.  (M?P  cross  sections). 
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Fig.  17.  Comparison  of  temperature  and  density  profiles  at  10  nsec 
including  radiation  transport  with  temperature  and  density 
profiles  omitting  this  effect.  (MPP  cross  sections). 
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Fig.  18.  Temperature  and  density  history  of  the  rear  surface  of  the  plasma 
with  and  without  radiation  transport.  (MPP  cross  sections). 
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Fig.  20.  Energy  history  of  the  beam  when  radiation  transport  is  omitted 
(same  labels  as  Fig.  15). 
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